H75-T5Cn^ 


'"(NaSfl^CR-120 520) STUDY OF KOLTIPLE HOLGGBAa 
FECOfiDING IH LITHIUM HIOEATB Final fieport, 

1 Oct* 1 97 3 - 31 Oct. 1974 (Georgia Inst, of 
Tech.) 55 p HC $4.25 CSCL 1 4E Unclas I 

V ^ G3/35 06364 J 

' ' ■ ' 1) 

Contract NAS8-30246 


STUDY OF MULTIPLE HOLOGRAM RECORDING IN LITHIUM NIOBATE 


T. K. Gaylord and W.R. Callen 
School of Electrical Engineering 
Georgia Institute of Technology 
Atlanta, Georgia 30332 


October 1974 


FINAL REPORT FOR PERIOD 1 OCTOBER 1973- 31 OCTOBER 1974 


Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 


US Department of Commerce 
Springfield. VA. 22151 


PRICES SUUEa TO QUiiGE 



Performed for 

National Aeronautics and Space Administration 
George C. Marshall Space Flight Center '' 
Marshall Space Flight Center, Alabama 35812 



SCHOOL OF ELECTRICAL ENGINEERING 
Georgia Institute of Technology 
Atlanta, Georgia 30332 


FINAL REPORT 
PROJECT NO, E-21-633 

STUBY OF MULTIPLE HOLOGRAM RECORDING IN LITHIUM NIOBATE 


by 

T. K. Gaylord and W. R. Callen 


RESEARCH CONTRACT NASS- 3 0246 
1 October 1973 to 31 October 1974 


Performed for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
George C. Marshall Space Fli^t Center 
Marshall Space Flight Center, Alabama 35812 



ABSTRACT 


The results of detailed experimental and theoretical considerations 
relating to multiple hologram recording in lithium niobate is reported. 
The following problem areas are identified and discussed; 1) the angular 
selectivity of the stored holograms , 2) interference effects due to the 
crystal surfaces, 3) beam divergence effects, 4) material recording sen- 
sitivity, and 5) scattered light from material Inhomogeneitles . 


i 



TABLE OF CONTENTS 


I. 

II. 

III. 

IV. 

V. 

VI. 

VII. 
VIII 

IX. 


INTRODUCTION 

ANGULAR SELECTIVITY 
"Angular Selectivity of Lithium 
Niobate Volume Holograms" 

SURFACE INTERFERENCE EFFECTS 

BEAM DIVERGENCE EFFECTS 

RECORDING SENSITIVITY i 

"Volume Holographic Recording and 

Storage in Fe-Doped LiNbO^ using Optical Pulses" 

SCATTERED LIGHT EFFECTS 

"Light Scattering Induced Holograms 

in Lithivnn Niobate" 

MULTIPLE HOLOGRAM STORAGE 

RECORDED HOLOGRAM ANALYSIS 
"Calculation of Arbitrary-Order Diffraction 
Efficiencies of Thick Gratings with Arbitrary 
Grating Shape" 

SYSTEMS CONSIDERATIONS 

"Optical Msmories: Filling the Storage Gap" 



REFERENCES 



LIST OF FIGURES 


Figure 

1. Experimental Configuration for Measuring the Optical 
Holographic Storage Properties of Ferroelectric Crystals 

2. A Comparison of Angular Selectivity Experimental Data 
with Results from the Coupled Wave Theory 

3. Amplitude Range of Theoretical Angular Selectivity Due 
to Minute Thickness Variations 

4. Transmittance Factor, Tj as A Function of Crystal Thickness 

5. Effect of the Transmittance Factor on Angular Selectivity 

6. Required Writing Energy Density for Various Optical 
Recording Materials 

7. Comparison of Coupled-Wave and Matrix Theories 

8. Schematic of a Read-Write-Erase Optical Holographic 
Computer Memory 



I., INTRODUCTION 


A program to study high capacity recording in electro-optic crystals 
was undertaken. This study focused on the problems associated with very 
high capacity storage (through multiple hologram superposition) and the 
use of lithium niobate as the recording mediimi, 

A nvimber of important problem areas were identified and studied. 

These included; 1) the angular selectivity of the stored holograms, 2) 
interference effects due to the crystal surfaces, 3) beam divergence effects, 
4) material recording sensitivity, and 5) scattered light from material 
inhomogeneities. 

Single hologram and multiple hologram recording experiments were per- 
formed on a vibration isolation table with an argon laser using the experi- 
mental configuration shown in Fig. 1. Holograms were analyzed experimentally 
by reading them with laser light of 488. Onm, 514. 5nm, and 632. 8nm wavelength. 

To properly interpret the above experimental results, a number of 
analytical studies were initiated to provide predictions of the read-out 
parameters of the volume holograms. 
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FIGURE 1 


EXPERIMENTAL CONFIGURATION FOR MEASURING THE OPTICAL 
HOLOGRAPHIC STORAGE PROPERTIES OF FERROELECTRIC CRYSTALS. 




II. ANGtFLAR SELECTIVITY 


Because of the three-dimensional nature of the crystal, holograms 
recorded in lithium nlobate are volume holograms. These holograms are 
produced directly (without processing) by the interference of laser beams 
of the appropriate wavelength intersecting in the crystal. The volume 
nature of this holographic storage is especially interesting since it 
indicates the possibility of very high capacity information storage. 

Volume (thick) holograms exhibit a number of properties in addition 
to those possessed by two-dimensional (thin) holograms. Among these pro- 
perties is angular selectivity—the need for the reference beam to illu- 
minate the hologram at a precise angle in order to achieve reconstruction. 
Illumination outside of this angular corridor produces a rapidly decreasing 
intensity of the reconstructed data. To perform these diagnostic experi- 
ments plane wave holograms (caused by the interference of two plane waves) 
were used. Even though this is a special case, a general hologram may be 
constructed by the superposition of an infinite number of plane wave holograms. 
Reconstruction of a voltane hologram Cl] is possible (with half maximum 
diffracted power or greater) only over the range of wavelengths given by 

AX L 

— ^ cot e - (1) 

X 

and only over the range of angles given by 

A6 ~ A(i^^)n(x,P) ^ , 

where and ^9 are centered about the writing wavelength 
angle, 0. The thickness of the hologram is given by t, L is 


( 2 ) 

and the writing 
the fringe 
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spacing of the fundamental grating (L = ^/Zsing), n(x,P) is the appropriate 
index of refraction for the probing beam wavelength and polarization, and 
A(t^) is the angular selectivity coefficient (approximately equal to unity). 
These types of properties are in actuality just manifestations of the 
increased storage capacity of the volume storage medium. 

The experimental configuration shown in Figure 1 was used to measure 
the angular selectivity and compare it with the theoretical value. Thicknesses 
from 1mm to 5mm were tested. Both iron-doped and nominally pure crystals 
of lithium niobate were used. The half-power angular widths were found to 
be in agreement with theoretically predicted values. These results show 
that the theoretical angular packing density of multiple holograms is 
achievable! This important finding was reported in Appl, Phys. Letters [2]. 
This publication gives experimental details as well as results and is repro- 
duced here for completeness. 

Additional experimental results have been obtained since Ref. 2 was 
published. These show a much more detailed comparison of theory and experi- 
ment for a broader range of reading angles. These data are shown in Pigs. 2 
and 3. 
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FIGURE 3 , AMPLITUDE RANGE OF THEORETICAL ANGULAR SELECTIVITY DUE TO MINUTE THICKNESS 
VARIATIONS. The thickness used in producing the plain curve is that 
needed to obtain a minimum value of the transmittance factor on the Bragg 
angle (t=0,3). The crossed curve corresponds to maximum transmittance 
(t’=1*8). The experimental data is seen to occur within these extremes. 


III. SURPACE INTERFERENCE EFFECTS 


Many analyses of hologram reconstruction do not account for boundary 
(surface) reflections. The diffraction efficiency results may be corrected 
to include boundary reflections by multiplying by a transmittance factor. 
This factor, T, is developed by us in Ref. 3, a copy of which is included 
in this report. This factor is the same as the transmittance factor derived 
by Kogelnik and given as Eq. (8) in Ref. 4, but with \id in that equation 
replaced by the argument of the sine function in Eq. (27) in our paper [3], 
From these results, we find that boundary reflections produced by the 
surfaces can considerably change the diffraction efficiency. The change 
can be an Increase or a decrease depending on whether the transmittance 
factor is greater or less than unity. This effect has been studied by 
Cohen and Gordon [5], For the grating parameters used here, t is typically 
in the range 0.70 to 1.20. In practice, the boundary reflections can be 
eliminated by antireflection coatings on the surfaces of the gratings. 

Calculated results showing the thickness dependence of t for lithium 
niobate are shown in Fig. 4. The effect of t as a function of reading 
angle is shown in Fig. 5. 
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FIGURE 5. EFFECT OF THE TRANSMITTANCE FACTOR ON ANGULAR SELECTIVITY 




IV. BEAM DIVERGENCE EFFECTS 


The effect of writing beam divergence was studied. First, a hologram 
was written with a frequency doubled Nd:YAG laser = 530nm) with a beam 
divergence of 0,143°. Second, a hologram was written with an argon laser 
(X = 514,5ran) with a beam divergence of 0-035°. The angular selectivities 
of these holograms were then measured with a He~Ne laser. The half-power 
angular width in both cases was found to be within 10% of the theoretical 
value derived in our paper [2']. Thus, it was concluded that writing beam 
divergence has little effect on the final read-out process. 




V. RECORDING SENSITIVITY 


Recording materials must possess a number of important characteristics 
to achieve the high storage capacities that have been predicted for optical 
memories. These requirements on the optical recording material include: 

1. High sensitivity— It is desirable that only a small amount of optical 
energy per unit area be needed to record the hologram of a data page« Table 
3 in Ref. 6 (reproduced later in this report) lists the necessary writing 
energy densities for a number of recording materials. For a practical system 
an energy density of about 1 millijoule/cm or less will be needed. 2. Large 

diffraction efficiency— Diffraction efficiency is the fraction of the reading 

\ 

light (reference beam) that is diffracted into the reconstructed data beam. 

It must be possible to record a single hologram with a large diffraction 
efficiency, so that in practice many holograms may be recorded at a single 
location, each with an equal share of the total maximum diffraction efficiency. 
Therefore, it is desirable to have the maximum diffraction efficiency as 
close to 1007. as possible. 3. Erasable and rewritable— For a rapid cycle 
read-write-erase memory system, it must be possible to continuously alter 
the stored data in the memory without encountering any degradation in the 
material characteristics. 4. Long lifetime of stored information— Stored 
data should persist for long periods of time before having to be refreshed. 
Ideally, storage should be permanent. 5. Non-volatile storage— Data should 
remain recorded in the memory in the absence of system power. 6. Nondestruc- 
tive readout— It should be possible to perform an essentially unlimited 
number of read operations without degrading or altering the stored data. 

7. Three dimensional storage — To achieve very high capacity storage, the 
information should be stored in thick (volume) holograms. Together with the 
requirement of high diffraction efficiency, this means that the hologram 


should be a thick phase (nonabsorbing) hologram. 8. High resolution-^ 

The storage material obviously must be capable of recording the very fine 
(wavelength size) variations of the Interference pattern produced by the 
intersection of the object and reference beams. 

Considering all of the above material requirements, the phot ©refractive 
materials (optically induced changes in index of refraction) appear to be 
especially promising. These materials, often ferroelectric crystals such 
as lithium niobate and strontium barium niobate (SEN), have been consider- 
ably developed and improved. For example,, in the first use of lithium 

niobate as a recording material in 1968 a writing energy density of approxl- 

2 

mately 100 joules/cm was required ^7]. Less than six years later, doped 
versions of lithium niobate have now been shown in this work to exhibit 

energy densities of 2 millijoules/cm ! We have announced this im- 
provement in sensitivity of almost 5 orders of magnitude in Applied Physics 
Letters [8]. This article is reproduced here for completeness. Fig. 6 
depicts this recent jump in sensitivity with respect to other potential 
recording materials. In addition, recent work by von der Linde et al. [9^] 
indicates that even higher sensitivities are possible in lithium niobate! 




VOLUME holographic RECORDING AND STORAGE IN Fe-DOPED LiNbO 3 USING OPTICAL PULSES 


PRADEEP SHAH, T A RABSON AND F K TITTEL 

Deoartment of Electric Engineering Rice Univ. Houston Texes 
T K GAYLORD 


paper INTENTIONALLY OMITTED 






Writing 

Energy Density 


Material 

Type of Material 

( |oules/em^ } 


°2D 

Ferroelectrie-Photoconductive 

1 X 10"® 


Malachite Green: 


2 X 10"5 


Sucrose Benzoate 

Thermoplastic 


Agfa 8E70 

Photographic 

2 X 10“® 


Kodak 649F 

Photographic 

7 X 10“® 



Ferroelectric-Photoconihjctive 

1 X 1(t3 


^LiNbO^:Fe 

Photorefractive 

2 X 10-3 

/ 

^'a75®^0.25''^*'2°6 

Oichromated 

Photorafractive 

6 X 10t3 


Gelatin 

Photochemical 

9 X 10“^ 

/ 

Ca F 2 : Ce 

Photochromic 

lx 19“^ 

/ 

KCL:Na 

Photochromic 

1 X 10“^ 

1 

Rice Univ. 

Mn Bi 

Magnetooptic 

3 X tO-2 

Georgia Tech 

TejjGe^ASg 

Amorphous Semiconductor 

5 X 10 ^ 

1974 

Gd 1 G 

Magnetooptic 

9 X 10“2 

1 

Eu Q 

Magnetooptic 

9 X 10~2 

1 

Na F 

Photochromic 

9 X 10“2 

1 

Co-P-Ni-Fe 

Magnetooptic 

1 X 10~^ 

\ 

SrTiO^tNhMo 

Photochromic 

2 X 10“^ 

\ 

Ba Ti Oj 

Photorefractive 

2 X 10“^ 

\ 

Mn Al Ge 

Magnetooptic 

3 X iir^ 



Amorphous Semiconductor 

5 X 10“^ 

8 X 10“^ 


^LiNb 03 :Fe 

Photorefractive 


KBr 

Photochromic 

1 


Cu^Hg 

Thermoplastic 

3 


BaNaNbgO^g 

Photorefractive 

5 


«Wl2 

Photorefractive 

10 


^‘'0.75®®a25*^‘’2°6 

Photorefractive 

14 


LiMbO^ 

Photorefractive 

100 
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VI. SCATTERED LIGHT EFFECTS 


Scattered light during hologram reconstruction has been recognized as 
a problem for high capacity storage in lithium niobate flO], 

We have reported [11] the presence of cones of diffracted light upon 
illumination of previously laser-exposed crystals of lithium niobate. These 
diffraction cones are. shown to result from the Internally recorded inter- 
ference pattern (hologram) resulting from the interference of the original 
incident laser beam with light scattered from material inhomogeneities, 
diffraction cones are observed in iron-doped lithium niobate crystals that 
were exposed to a single laser beam and in crystals that were exposed to 
two superposed laser beams (i.e., during conventional holographic recording). 
In the two beam case, the diffraction cones are present in addition to the 
first order diffracted beam when the conventional two beam thick hologram 
is reconstructed. The diffracted cones produce the impression of scattered 
light during hologram reconstruction, an effect that has previously been 
reported in transition metal doped lithium niobate [10]. 

The diffraction cones, which have their apex in the exposed region of 
the crystal, are observed as rings (referred to as "scattering" rings, or 
diffraction rings) when a screen or a piece of film intersects the cone of 
light. Figure 1 in Reference 11 shows two typical diffraction ring patterns. 
For the single beam case, the observed results in lithium niobate are effec- 
tively the same as the experimental observations of Moran and Kaminow [12] 
for polymethyl methacrylate (PMMA) , which had been exposed to ultraviolet 
laser light. 

The presence of diffracted cones of light represents a possible limi- 
tation of heavily iron doped lithium niobate for data storage applications 



because optical power is lost into the scattering induced diffraction cones 
that could otherwise be used to increase the diffraction efficiency and thus 

the total bit capacity of the two beam grating hologram. However, it has 
already been shown by Phillips, Amodei, and Staebler [101| that the "scattered" 
light may be erased 1) by illumination with uniform incoherent light or 2) by 
writing additional superposed holograms at new angles. In the latter case, 
"scattered" light from the previous holograms tends to be erased. 

Our Ref. 11 is reproduced here for completeness. 
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VII. MULTIPLE HOLOGRAM STORAGE 


In the last four sections, problems associated with multiple hologram 

storage have been discussed. The experimental results presented have been 

obtained using the basic experimental configuration shown in Fig. 1. Basic 

diagnostic experiments were performed by storing both single holograms and 

by storing multiple holograms [13^ at a single location. The theoretical 

8 2 

storage density of two dimensional (thin) holograms is 4 x 10 bits/cm (one 

bit per square area one wavelength on a side) whereas in three dimensional 

12 1 

volume (thick) holograms the theoretical storage density is 8 x 10 bits/cm' 
(one bit per cube volume wavelength on a side) [14]. Obviously for truly 
high capacity storage, thick holograms (such as in optical crystals) need to 
be used instead of thin holograms (such as in photographic emulsions or 
metal films). Holographic memory systems have been described that utilize 
three-dimensional storage [15 3- These systems superpose many holograms at 
a single location inside the thick recording medium by using a different 
reference beam angle for each hologram. The superposition of multiple 
holograms at a single volume location introduces the additional problem of 
writing new holograms in that volume without affecting those already there. 
When lithium niobate is used as the three dimensional storage material, this 
problem may be solved by the application of an external electric field [I 63 , 
[I 73 . This greatly increases the sensitivity for writing while the sensi- 
tivity for erasure remains unchanged at a much lower value. Thus, as a new 
hologram is written, the other holograms at that location are only slightly 
erased. Work is presently underway in our laboratory to duplicate these 
electric field effects. 



VIII. SECORDED HOLOGRAM MALYSIS 


A method for analyzing the diffraction efficiency of thick, lossless 
transmission holograms in lithium niobate was developed. In lithium niobate 
and similar ferroelectrics , the literature assumes the induced changes in 
index of refraction are sinusoidal in nature, like the two beam plane wave 
interference pattern. The diffraction efficiency can be predicted for the 
sinusoidal case [l^* In actual fact, the index of refraction variation is 
probably not sinusoidal due to the obviously nonlinear writing characteristic 
(diffraction efficiency versus exposure), which is experimentally observed. 

We have developed [S'] a method for calculating arbitrary-order diffrac- 
tion efficiencies of thick, lossless transmission gratings with arbitrary 
periodic grating shapes. For illustration, numerical values of the dif- 
fraction efficiencies at the first three Bragg angles were calculated for 
sinusoidal, square wave, triangular, and sawtooth gratings. The complete 
details of this method are expounded in Ref. 3, which is duplicated in this 
report. Also a comparison of our method to an extension of the Burckhardt 
matrix method [18] is presented in Fig. 7. Our method was determined to be 
20 times faster on the computer! 
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A method for calculating arbitrary-order diffraction efficiencies of 
thick, lossless transmission gratings with arbitrary periodic grating 
shapes has been developed. This represents an extension of previous 
work to nonsinusoidal gratings and to hi^er-order Bragg angles. A Fourier- 
series representation of the grating is employed, along with a coupled- 
mode theory of diffraction. For illustration, numerical values of the 
diffraction efficiencies at the first three Bragg angles are calculated 
for sinusoidal, square-wave, triangular, and saw-tooth gratings. Numerical 
results for the same grating shapes with the same parameters are also 
calculated for comparison, by extending Burckhardt's numerical method for 
analyzing thick sinusoidal gratings . The comparison shows that the coupled- 
mode theory provides results with relative computational ease and results 
that are in agreement with calculations obtained by extending the more -rigorous 
Burckhardt theory to nonsinusoidal grating shapes and to higher-order Bragg 
angles . 
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It is veil knowit that thick dielectric diffraction gratings differ 

from thin gratings in a number of important ways . Among these are the 

capability of high diffraction efficiency,^ wavelength selectivity,^ 

1 2 

angular selectivity, and reduced noise. These give rise to the use of 

thick gratings as highly efficient diffraction gratings, narrow-band 

3 4 

spectral filters, thick-grating optical components, such as lenses, and 

. 2 
imaging systems capable of spectral resolution of extended objects. In 

the field of integrated optics, thick gratings may be used as diffraction 

5 

gratings for surface guiding of waves, for thin- film distributed-feedback 

lasers,^ for frequency-selective grating reflectors in thin-film lasers,^ 

for grating couplers for launching single-mode ligjit waves into thin-film 

waveguides, and for electro-optic grating deflectors and modulators . 

In addition, thick (volume) holograms may be regarded as recordings 

of an infinite number of thick gratings. Thick holograms have attracted a 

11 

great deal of interest by their use in high-capacity information storage, 

12 13 

in color holography, and in white-light reconstruction of holograms. 

The diffraction of a plane wave by a thick sinusoidal grating at or 

. 14 1 

near Bragg incidence has been considered by Burckhardt and by Kogelnik. 

Burckhardt has treated this case by solving the exact electromagnetic 

boundary -value problem and has obtained numerical results with a digital 

computer to determine the eigenvalues of a matrix and to solve the resulting 

set of linear algebraic equations. Kogelnik has obtained a closed-form 

expression for the diffraction efficiency at the first-order Bragg angle, 

by employing a coupled-wave theory. Coupled -wave theories also have been 

1 S 

used successfully in the treatment of light diffraction by acoustic waves. ’ 
Recently, Chu and Tamir^^ treated this problem by using a guided-wave 



technique. They assumed sinusoidal modulation of the relative permittivity 
by the sound wave. Their treatment was based on a rigorous modal approach, 
utilizing the interrelationships between the characteristic mode and the 
coupled-mode representations. With their method, not only the diffraction 
efficiency at the first order but also that of any higher order can be 
obtained . 

In this paper, Chu and Tamir's approach is extended to examine the first 

and higher-order diffraction efficiencies of thick, arbitrary- shape gratings. 

Because of the periodicity of the grating, a Fourier- series representation 

of the grating is employed. The gratings are assumed to be lossless. The 

reflections at surfaces of the gratings are at first neglected in the 

derivation, because, in practice, these can be eliminated by antireflection 

coatings. When surface and internal reflections are present, the results 

18 

are corrected by a multiplicative transmittance factor. For illustration, 
numerical values of the diffraction efficiencies at the first three Bragg 
angles are calculated for sinusoidal, square-wave, triangular, and saw-tooth 
gratings. For comparison, numerical results for the same grating shapes with 
the same parameters are also calculated by extending Burckhardt's numerical 
method. The comparison shows that the results from these two methods are in 
close agreement and that the present method is computationally simpler and 


more efficient . 



THEORETICAL ANALYSIS 


The model for a thick periodic grating can be described by Fig. 1. The 
X axis is chosen in the plane of incidence and parallel to the surfaces of 
the medium, the z axis is perpendicular to the surfaces of the medium, and 
the y axis is perpendicular to the page. For convenience, the fringe planes 
of the grating are assumed to be perpendicular to the Surfaces of the medium 
and to the plane of incidence. The grating vector K is, therefore, parallel 
to the X axis. Thus, for lossless periodic gratings, the fringes of the 
grating can be represented by a spatial modulation of the relative dielectric 
constant, 

CQ 

Cj-(x) = S (x)/®o - Cro "*■ ^ [e^j^cos(hKx) + e^j^sinChKx)] , (1) 

h-1 

where K = 2 n/L, L is the period of the grating, is the average value of 
e^, and and are the spatial modulations of e^, the subscripts c, s 
and h denoting the quantities connected with the cosine gratings, the sine 
gratings and the h th-harmonic grating, respectively. Corresponding to the 
distribution of the relative dielectric constant, the distribution of the 
refractive index of the medium is 

O) 

n(x) = + 2 [n^^cos(hKx) + in (hKx) ] , (2) 

h=l 

/ 

where n is the average refractive index of the medium, and n , and n , are 
o ch . sh 

the spatial modulations of n. 

The electric field of the incident wave is assumed to be polarized 
perpendicular to the plane of Incidence (H mode) and is of the form 



exp[ j(Tj^x + • The wave propagation in the grating can be described 

by the scalar wave equation 

\ 

[ 7 ^ + k^e^(x)] E(x,z) = 0, (3) 

where k = 2 rr/ \ \ is the £ree-space wavelength of the incident plane wave, 
and E(x,z) is the complex amplitude of the y component of the electric field, 
which is independent of y. 

14 

Eq. (3) has been solved by Burckhardt, using separation of variables, 

an infinite-series solution for the x-dependent equation, and a matrix method 

to solve the eigenvalue problem associated with a truncated set of the resulting 

infinite system of equations . This approach has been used to obtain numerical 

19 

results for sinusoidal gratings. Kaspar has extended Burckhardt 's method to 

find the diffraction efficiency for nonsinusoidal absorption gratings. He has 

pointed out that when the absorption is strong, the phase-grating contribution 

to the diffraction efficiency is very small. 

17 

Chu and Tamir have shown that the field inside the grating can be 
described in terms of coupled modes if the modulations of the relative dielectric 
constant are very small. In general, in addition to the zeroth mode, many 
higher-order modes are excited, because of the presence of the grating. For 
an incident wave (zero-order mode) of wavelength X , and at an angle 9^ , the 
fundamental grating will diffract this wave if the Bragg condition, m X=2L sin 9 q , 
is satisfied or nearly satisfied. For this particular wavelength and angle, the 
harmonic gratings may or may not produce diffraction, depending on whether or 
not their corresponding Bragg conditions are satisfied or nearly satisfied. 

These diffracted modes of the fundamental and the harmonic gratings propagate 
in the same direction, 

( 



The dimensionless quantities 


Ich ■ 2< 

L 

h X 

) ^ 

^ ech’ 

h = 1, 2, 3, . . . 

(4) 

-sh = 

L 
h \ 

^ *sh ’ 

h = 1, 2, 3, . . . 

(5) 


are c 


ailed the effective-modulation indices Because and would 


’sh 


typically be lO”^ or smaller, and q^^ are small even if L is many times 

as large as For example, if 0^= 5° and = 10 then q^^ = 0.0066. 
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When q is very small compared to unity, it can be shown ’ that two 
coupled-wave equations and therefore two modes are sufficient to describe 
the coupling effects when the incident angle is equal to or near the Bragg 
angle. Therefore, for an incident wave of wavelength \ and at an angle 8^ , 
the electric field inside the grating can be written as the sum of the 
fundamental mode and an arbitrary mode 


E(x,z) = S^(z) exp(jT^x) + S^(z) exp(jT^x) 


< 6 ) 


where X| ^ zeroth-mode and the m th-mode (with respect to the 

fundamental grating) transverse wave numbers, respectively. The continuity 
of the electric field at z = 0 and the Floquet theorem require that Tl^ - 1)^ 

- k sin 6^ and = 11^ *• 2mrr /L. The tilde ~ will henceforth be used to 

denote the quantities in the dielectric medium when the gratings are present. 
The integer subscript m represents the m th-order diffraction. The integer 
h represents the h th-harmonic grating. Diffraction occurs when the h th- 
harmonic grating satisfies or nearly satisfies the Bragg condition A. = 
2(L/h) sin9^, where rtij^ represents the m^^ th mode (with respect to the h th- 
harmonic grating) excited due to the h th-harmonic grating. Exact-Bragg 
conditions occur when m^ is equal to m/h where h divides evenly into m. 
Near-Bragg conditions occur for the wavelength X 1) when the angle of 

/ 



incidence is near, but riot equal to 6^ , and/or 2) when the value of m is 

large, and h divides nearly evenly into m, so that m/h is almost an integer. 

Thus, in Eq. (6) represents the total amplitude, together with the 

propagation factor in the z direction, of the diffracted mode due to all of 

the gratings that satisfy or nearly satisfy the foregoing Bragg condition. 

At the boundary z = d, propagates at the angle 0^ , whereas propagates 

at an angle 9 , which is determined by 
m 

9^ = - sin”^(-^) = - sin ^(sind^ - (7) 

The diffracted modes due to the gratings that are far from obeying the Bragg 

condition are assumed to be negligibly small compared with and 

Therefore, the interaction between ^ and *S can be characterized by the 

o m 
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coup led -mode equations 
dS _ m 

T-^ - j| ”§ - J[2^(C . + jC . )1S = 0, (8) 

dz *^^o o chmj^ shm^^ ra ^ ^ 

dS _ ^ m 

S - j[2(C . ' jC . )1‘§ = 0, (9) 

dz '’^m m ® ^ ' 

where § and ^ are the longitudinal wave numbers inside the medium when the 
o m 

— 1 - — 2 

gratings are absent. They are given by k(c^^)^cos cp and {k [k 

^ 2 % 

(e^^)^sin cp - (2mn/L)] ] > where cp, the refraction angle in the medium, is 
-1 h- 

given by q) = sin [(sin ® *] • The bar notation will henceforth be 

used to denote the quantities Inside the medium when the gratings are absent. 

For a given value of the integer m, the subscript h may be any integer that 
divides evenly or nearly evenly into m, provided the corresponding h th-harmonic 



m 


grating exists. The symbol S denotes the summation over all of these possible 

^ 17 21 

values of h. The coupling coefficients in Eqs. (8) and (9) are given by ’ 


C , 2! 

ch 


[■ 


^o (2) (n^-1)! 


.l2, JL_ s2, A 
^ A/h ^ ^‘*ch^ ’ 


( 10 ) 


c , - [ 

“■’"h r + 5 


o' ^m (2) (m^-1)! 


.i2/ JL_ ^2. A 

^ ^ L/h ^ * 


( 11 ) 


If the grating does not exist, there is no coupling between 

S and S and therefore no diffraction. Under this condition, only Eq. (8) 
has physical significance. It represents the propagation of the fundamental 
mode (incident wave) inside the medium. 

The solutions of Eqs* (8) and (9) are of the form 


\(z) = A^exp(jl^z) + B^exp(j§^z), 


(12) 


= A^exp(J5^z) + B^cxp(j§^z). 


m 


m 




(13) 


The wave numbers | and § can be found directly by substituting Eqs. (12) and 


m 


(13) into Eqs. (8) and (9). They are 


§ 


^o”^ +^.^o ^mv2 


o,m 


[( 


^ 9 ^ 9 h 

)"* + ( 2 C . )^ + ( S C . )^]^, 

h '"V h ^ 


(14) 


where the + sign corresponds to and the - sign to Tbe constants A^, 

B , A , and B are determined by the boundary conditions and 
o m m ■' 


ra 


[E (C , - jC . )]A - ( § - e )A = 0, 

shnij^ 0 ^o ’tn m 


( 15 ) 





(16) 


m 


)B„ = 0. 


h 


nij^ -’'^shni^^'-' o ' ^in 'm 


which are obtained from Eqs , (9), (12), and (13). To specify the boundary 
conditions, the amplitude of the incident wave is assumed to be unity at z=0 
so that, from Eq. (12), 


S (0) = A + B = 1, 
o o o 


(17) 


Initially, the amplitude of the diffracted wave is zero. Therefore, evaluating 
Eq. (13) at z = 0 gives 


S (0) = A + B = 0. 
m m m 


(18) 


Solving Eqs, (15), (16), (17), and (18) for A^, B^, A^, and B^ gives 


A = 


^Q-^m 
?o ■ ^m 


- — — — 0 9 9 

1 <Sq - sj t(l„ - 4 [( H- ( 

^ [<C -|')^+4[(EC )^ + ( EC )^])^ 


m 


(19) 


B = 




”ra 


m 


- 5o) 


§„)^ + 4 [ ( EC 


m 


chmh 


) + ( SC 


h 




[(§ -5)^+4 [(EC )^ + ( EC )^33^ 


( 20 ) 


A = 
m 


- B = 
m 


lu m 

( SC . ) - j ( SC . 

h h 


— — 9 ^ 9 ^ 9ir 


( 21 ) 
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For the exact-Bragg condition, and 0^ “ ®n ‘ Hence, Eqs, (lA), (19), 


-m 


(20), and (21) become 


I =g t[(2C, ) + ( 2 C , 


(2Z) 


A = 


A = 
m 


1/2 = 


- B - ^ 
m 2 


m 

(SC. 
h <=’''% 


m 


) - j( 2 C . ) 

h 


m 


m 


C( 2 C )" + ( 2 


h 


shmh 


(23) 


(24) 


Thus, the transmitted and the diffracted modes are 


m 


m 


S^(z) = exp(j?^z) cos{[ ( 2 + ( 2 




(25) 


- ^ m 9 k 

S(z) = j2A exp(jf z) sinCC ( 2 C . ^ )^ + ( S C , ) V , (26) 


m 


m 


“•“Y' 


shnlh 


where and are given by Eqs . (10) and (11) with and 

^ ”\i "*h 

is given by Eq. (24). Eq. (26) is the general formula for the m th-dif fracted 
mode due to any periodic grating when the incident angle of the zeroth mode 
satisfies the Bragg condition m \= 2L sin 0^ . The diffraction efficiency for 
the m th order of diffraction is defined as 


. Sjd) S*(d) 

DE ^-!!! 1— ■ 

s^(0) s^(0) 

and thus for exact-Bragg conditions 


(27) 


m 


, m ui - , 

DE = sin {[(2C, ) + (2 C , ) 1 d} , 

h h 


m 


(28) 


where the asterisk * denotes complex conjugate. Upon substituting Eqs. (10) 

3 / .) 



and (11), with into Eq. (28) and performing some algebraic manipulations, 

we find that 


DE = sin 
m 


[(f 


[• 


(2) (mj^-1)! (h) 






(2m,-l) 






m 

+ C S 


[• 


ra, 

^ (2) “ (V^>- 


3^2<v1) 




(2m^^-l) 


X - (e^o)coscp 


, hi 

y\. 


(29) 


Eq. (29) is the general expression for the diffraction efficiency at the m th- 
order Bragg angle for a periodic grating of arbitrary grating shape. For example, 
the first-, second-, and third-order diffraction efficiencies for a grating, whose 
dielectric constant profile can be expressed as a Fourier sine series, are 


DE, 


= sin^- 




2 \ (e^^) cos «p 




t 2 / ^2 

^ (e.i) 


,, . 2rr "'•sr , T nd 

DE, = sin [[- -j — + e ,] I 

2 2 r 2 Ue r COS9 


}. 


(30) 


(31) 


and 


®sl^ 

DE, = sin tC 7~~ 6 

^ 16 V 


rrd 


s3^ 


2 cos<p 


}■ 


(32) 


In Eqs. (30), (31), and (32), only the Fourier grating components 6^2* 

0 are required to evaluate the diffraction efficiencies DE , DE«> DE , 

s3 l z o 

Table I gives these Fourier components, normalized to the amplitude of the 
fundamental grating, for gratings having sinusoidal, square-wave, triangular, 

and saw-tooth dielectric constant profiles. Note that the sinusoidal, square- 



wave, and triangular grating shapes can each be represented by a Fourier cosine 
series also. In this case, the resultant diffraction efficiency expressions 

contain only e^ 2 ’ ®c3* ”ch ^sh "o’ true in 

most cases^, it can be shown that ® ^"o"ch ®sh ” ^"o"sh' Therefore, 

k 

with (s^q) = Eq- (29) becomes 


m 


DE = sin 
m 






2 L 


(m^-1)! (h) 


(VlV (2n^,-X) 


2(m-l> , ,”h 

^ n(n_) ( n^i,) j 


coscp 


m 

+ { 2 [ 

h 


-n3 


2 L 




sh^ 


(mj^-1)! (h) 


(\-iy 


coscp 


5% 

1 ; '"J. 


(33) 


The results calculated with Eqs. (30), (31), and (32) do not agree with 
those calculated by use of Burckhardt's matrix method. This is because Burckhardt 
takes the boundary reflections into account, whereas they are not included in 
the foregoing derivation. Our results may be corrected to include boundary 
reflections by multiplying the diffraction efficiency by the transmittance 
factor, 


(1-R) [1+2RCOS 


(2pd)+R^ 


T 

m 


(l-R^)^ + 4R^cos^(2v^|jd)-h:os^(2gd)] - 4R(l+R^)cos(2v^d)cos (2pd) 


(34) 


m 

C( s c 

h 


chmj^ 


y 


2 2 k 

where R = sin (6^-^tp)/sin (6^+<p), p = 2 n (e ^) (coscp )A, and y 

+(SC, ) i evaluated with t = § for exact-Bragg conditions. This factor 

^ shm^ > ^m ^o 

is the same as the transmittance factor derived by Kogelnik and given as Eq. (8) 


in Ref. 18, but with vd in that equation replaced by the argument of the sine 
function in Eq. (28) of this paper. This allows generalization to higher 
diffraction orders and nonsinusoldal gratings. 



RESULTS AND DISCUSSION 


, The coupled-wave analysis In the preceding section was numerically 
implemented on a UNIVAC 1108 computer and calculations were performed for 
gratings having sinusoidal, square-wave, triangular, and saw-tooth distributions 
of the dielectric constant. Table II gives numerical values f,or the diffraction 
efficiencies at the first-, second-, and third-order Bragg angles for these 
gratings. These results represent DE^ DE 2 and DE^ as obtained 

from Eqs. (30), (31), (32), and (34) with 2.3225 (value used in Refs. 14 

-4 

and 18) and e ,= 10 e • The fundamental spacing of these gratings is L = 
s X iro 

JL o 

3.630 p,m (resulting from recording with two beams of 632.8 nm at 0^=T5.O ), 

For comparison, the results obtained by extending Burckhardt’s numerical method 

(matrix method) to nonsinusoidal gratings are also shown in Table II. These 

results were Calculated by programming Burckhardt's method on a UNIVAC 1108 

computer and using the UNIVAC Math Pack subroutines to solve the eigenvalue 

problem and the set of linear algebraic equations. Table II shows that the 

results of these two methods are in close agreement; the deviation between 

these two methods does not exceed 2.8% for diffraction efficiencies larger 

than 5x10 ^%. Diffraction efficiencies smaller than 5x10 ^% are less significant 

physically because the corresponding low-level diffracted intensities are 

•8 

difficult to measure. Diffraction efficiencies of less than 5x10 % have been 

listed as zero in Table II. In addition to the results in Table II, we have 
performed calculations for other grating thicknesses (15)Xm, ISOOjXtn, and 

2000^n0 and other fundamental grating spacings (1.222p,m and 1.822|xm). We found 
that the deviation between the coupled-wave analysis and the matrix analysis 
does not exceed 6.7% for any case with a diffraction efficiency larger than 



5x10 Typically, the percentage deviation is a few tenths of one percent. 

Although Burckhardt's numerical approach is rigorous, a number of 
mathematical problems such as truncation of the matrix and discarding of large 
positive eigenvalues must be overcome. A discussion of these is included in 
Ref. 14. In addition, another mathematical difficulty associated with the 
Burckhardt method, encountered in the present work, is a singularity that 
arises in the process of solving a set of linear algebraic equations. For 
pure phase gratings, Eq. (9) in Ref, 14 is real and symmetric. When the incident 
wave is at the Bragg angle, pairs of equal elements are introduced on the 
principal diagonal of the matrix in that equation. Thus, when the modulation 
amplitude is small, pairs of equal eigenvalues are usually induced. This results 
in a singularity in the matrix in Eq. (34) in Ref. 14; therefore, the equation 
is nonsolvable. For the parameters in the particular examples of Ref. 14, this 
problem does not occur because the modulation amplitude is large (0,0035 e^^). 
However, the modulation amplitude may, in practice, be very small (of the order 
10 ^ or smaller) and the singularity problem must, therefore, be overcome. A 
Mt&y to avoid the singularity is by shifting the incident angle by a negligible 
amount away from the Bragg angle. Physically, because the shift is negligibly 
small (10 ^ degrees was used here) , the incident wave can still be regarded as 
being incident at the Bragg angle. In the present method, a closed-form 
expression for the diffraction efficiency is obtained, and no mathematical 
difficulties arise in the process of calculation. The computer time needed in 
the present method is only about 1/20 of that needed with the extended Burckhardt 
method to perform the same calculations . 

From the results, we found that boundary reflections produced by the surfaces 



can considerably change the diffraction efficiency . The change can be an 

increase or a decrease depending on \diether the transmittance factor is greater 

23 

or less than unity. This effect has been studied by Cohen and Gordon. For 

the grating parameters used here, T is typically in the range 0.70 to 1.20. In 

practice, the boundary reflections can be eliminated by antireflection coatings 

on the surfaces of the gratings. We also found that the diffraction efficiency 

of a given higher order is mainly contributed by the corresponding higher- order 

Fourier component of the grating. The difference between the diffraction 

efficiencies for sinusoidal and nonsinusoidal gratings (having the same average 

and fundamental grating amplitudes) appears only in the higher- order diffractions. 

The higher-order diffraction efficiencies, however, very strongly depend on the 

grating shape. Also, for small grating modulations, the diffraction efficiencies 

at any order are very dependent on grating thickness ; they increase with 

increasing thickness. Marcuse^^ has suggested that, for small-amplitude thick 

nonsinusoidal phase gratings, the higher-order diffraction efficiencies might 

be estimated from the relative amplitudes of the spatial harmonics, consistent 

with the assumption of perturbation theory that only one Fourier component can 

satisfy the Bragg condition for a given wavelength incident wave. Our calculations 

show that this is true except when the amplitude of the harmonic grating (h “ m) 

is very small compared to the amplitude of the fundamental and the lower-order 

contributing harmonic gratings. In this case, the contributions from higher- 

order diffractions’ (h<m) are significant. In addition, we found that the agreement 

between the coupled-wave method and the matrix method is better when the h ^ ra 

25 

term is dominant over h<m terms. Rigrod has shown that for reflection gratings 
there is no correlation between higjier-order diffraction efficiencies and the 
corresponding harmonics of the index profile. The present results show that 



this is not true for transmission gratings. 

The present method can be used to analyze the diffraction efficiency of 

any thick periodic grating regardless of the dielectric constant profile (grating 

shape). The examples analyzed here have had even or odd symmetry. However, 

the method does not require any synmietry to exist, but only that the grating 

be periodic. From the gratings analyzed, different grating shapes have shown 

different distributions of higher-order diffraction efficiencies. This indicates 

the possibility that this type of analysis might be used in reverse to determine 

the grating shapes of thick hologram gratings such as those recorded in 

26 

ferroelectric crystals. Due to nonlinearities in these materials, a sinusoidal 

exposure does not necessarily produce a sinusoidal change in index of refraction. 

Depending on which of the possible physical mechanisms is operative in a given 

situation (such as drift of charge carriers or diffusion of carriers) different 

27 

grating shapes are generated. 

Further, the derivations in the preceding section have assumed that the 
grating medium is lossless, that the gratings are unslanted with respect to the 
grating boundaries (grating vector parallel to surfaces of medium), and 
that the incident wave is H mode polarized. If the medium is lossy, the results 
still apply except that the coupling coefficients are complex, and therefore 
the attenuation factors are implicitly contained in the expressions for the 
transmitted wave and the diffracted wave. The method presented, here can also 
be straightforwardly applied to the analysis of slanted gratings and to E mode 
polarization of the incident wave. 


jf ! 



CONCLUSIONS 

A simple method of calculating arbitrary- order diffraction efficiencies of 
thick transmission gratings with arbitrary periodic grating shapes has been 
presented. The analysis uses a coupled-mode theory to obtain a closed -fo2rm 
expression for the diffraction efficiency of an arbitrary order. This method 
provides results with relative computational ease and results that are in 
close agreement with those obtained by extending Burckhardt's numerical method. 
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FIGURE CAPTION 


Geometry of a thick grating with unslanted fringes. The 
spatial modulation of g is indicated by the line pattern. 




Fig. 1 


TABLE I. First three Fourier components for various relative dielectric 
constant profiles (grating shapes). Components are normalized to the 


amplitude of 

the fundamental 

grating, 




Grating 

Component 

Sinusoidal 

grating 

Square- wave 
grating 

Triangular 

grating 

Saw- tooth 
grating 


1 

1 

1 

1 

hi' "si 

0 

0 

0 

-1/2 

"s3^ "si 

0 

1/3 

-1/9 

1/3 


jb / 


TABLE II. Comparison of diffraction efficiency in percent at the first-, second-, and third-order Bragg angles 
for transmission gratings with boundary ref lections . and with the same average and fundamental Fourier grating 

components. The grating parameters are = 2.3225 (value used in Refs. 14 and 18), ®sl ” ^ro* 

“8 

L = 3.6303 p,m, and the wavelength X. = 0.6328jj|,m. Diffraction efficiencies of less than 5 x 10 % are listed as 

O.OOC-5), 








Grating 

thickness 

Diffraction 

order 


Diffraction Efficiency 
(in 7o with power of ten in parentheses) 


Sinusoidal 

grating Square-wave grating 

Triangular grating 

Saw-tooth grating 

(microns) 


Coupled-- 

wave 

Matrix Coupled Matrix 

wave 

Coupled- 

wave 

coupled- 

wave 



1 

1.64(-3) 

1.64(-3) 

1.64(-3) 

1.64(-3) 

1.64(-3) 

1.64(-3) 

1.64(-3) 

1.64(-3) 

10 

2 

0.00(-5) 

0.00(-5) 

0.00(-5) 

0.00(-5) 

0.00(-5) 

0.00(-5) 

3.81(-4) 

3.87(-4) 


3 

O.OOC-5) 

0.00(-5) 

1.14 (-4) 

1.14 (-4) 

1.26(-5) 

1.25(-5) 

l.l4(-4) 

L.13(-4) 


1 

1.57(-1) 

1.57(-1) 

l.57(-l) 

1.57(-1) 

1.57(-1) 

1.57(-1) 

1.57(-1) 

1.57(-1) 

100 

2 

0.15(-5) 

0.13(-5) 

0.15(-5) 

0.14(-5) 

0.15(-5) 

0.13(-5) 

2.49(-2) 

2.56(-2) 


3 

0.00(-5) 

0.00(-5) 

1.12(-2) 

1.12(-2) 

1.25(-3) 

1.23(-3) 

1.12(-2) 

1.12(-2) 


1 

1.2K+1) 

1.21(+1) 

1.21(+1) 

1.21(+1) 

1.21(+1) 

1.21(+1) 

1.21(+1) 

1.21(+1) 

1000 

2 

2. 53 (-4) 

2.46(-4) 

2.53(-4) 

2.46(-4) 

2. 53 (-4) 

2.46(-4) 

4.L6(+0) 

4.26(+0) 


3 

0.00(-5) 

0.00(-5) 

l,64(+0) 

l,64(+0) 

1.83(-1) 

l.S4(-l) 

1.64(+0) 

1.66(+0) 



IX. SYSTEMS CONSIDERATIONS 


A three-dimensional lithium niobate recording and storage system is 
shown schematically in Fig. 8. The systems aspects of such an optical 
recording scheme were thoroughly reviewed in this study. One of the results 
of this review was the publication of a state-of-the-art review [6"]. This 

article is reproduced in this report and is a self-contained review. An- 

/ 

Other result of this review was a change in our experimental reading system. 
An angular accessing system was developed and it is illustrated in Fig* 9. 
This system allows accurate and simple angular beam positioning without the 
need to rotate the crystal* 
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FIGURE 8. SCHEMATIC OF A READ-WRITE-ERASE OPTICAL HOLOGRAPHIC COMPUTER MEMORY. 
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FIGURE 9. ANGULAR ACCESS SCANNER FOR READING MULTIPLE HOLOGRAMS STORED AT A SINGLE LOCATION. 
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